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An analysis of low-energy charged pion-nucleon data from recent yr^p experiments is 
presented. From the scattering lengths and the GMO sum rule we find a value of the 
pion-nucleon coupling constant of = 0.0764 ± 0.0007. We also find, contrary to most 
previous analyses, that the scattering volumes for the P31 and P13 partial waves are equal, 
within errors, corresponding to a symmetry found in the hamiltonian of many theories. For 
the potential models used, the amplitudes are extrapolated into the subthreshold region to 
estimate the value of the E-term. Off-shell amplitudes are also provided. 
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I. INTRODUCTION 



The pion-nucleon interaction has been a fruitful source of knowledge of the strong interaction. The prop- 
erties of the baryon resonances produced in pion-nucleon collisions give strong support for the quark model. 
The subthreshold amplitude is related to the value of the pion-nucleon E-term which constrains models 
of nucleon structure. The pion-nucleon coupling constant provides fundamental input for the calculation 
of nuclear forces. Beyond these surface features lie more subtle indications of the nature of the strong 
interaction. 

We model the dynamics of this interaction by a coupled-channel Klein-Gordon equation whose potential 
is assumed to be the fourth component of a Lorentz vector. There are several advantages to this approach. 

First, unitarity, Coulomb corrections, multichannel effects, and hadronic mass splittings may be included 
in a natural way. 

Next, by solving the Klein-Gordon equation at the appropriate kinematical point, any observable can be 
calculated, even in the subthreshold region. This is especially easy in the case of the s-wave scattering with 
an exponential potential, for which an analytical result is available. 

By using such a model we are presented with an alternative approach to such quantities as the E-term. 
While the value of the S-term may well be more accurately determined by dispersion relations, by looking 
at it from the point of view of a potential theory the structure of the system is perhaps better revealed in 
the sense that the behavior of the subthreshold amplitude is directly related to the shape of the potential 
(and presumably to the distribution of the constituents of the pion and nucleon). 

A more subtle advantage is that models which take a broader view of the system, including higher energy 
data, have the problem that the model must be valid over the entire range. In this way the low-energy 
parameters have been largely determined by the data at high energies and the assumed dependence of the 
model for the low-energy behavior. Thus such features as the singularities in the scattering amplitude due 
to cuts coming from the range of the interaction depend almost entirely on the model assumptions. In 
the present technique, which does not claim to be a theory of the system, the low-energy regime can be 
investigated without recourse to data at higher energies. 

Finally, solutions in coordinate space allow us to develop an intuitive picture of the spatial structure of 
the interacting pion-nucleon system. 

Our Klein-Gordon model also has several drawbacks: The model is purely phenomenological at the 
hadronic level. Because it is a model based on static potentials, virtual particle production and annihi- 
lation and retardation effects are not explicitly included. As with all potential models, effects of crossing 
symmetry must be inserted. Because of the efficiency of the Jost calculation of the s-wave amplitude for the 
exponential potential (and the rapid cut-off properties of the Gaussian potential in r-space) we are able to 
incorporate (in a controlled approximation) the crossing symmetry condition that the isovector amplitude 
must vanish at the Cheng-Dashen (CD) point into the fitting procedure. The relativistic effect in the center- 
of-mass motion is taken into account only approximately. Because we describe only low-energy phenomena, 
we believe that these defects are outweighed by the model's strengths. 

In a previous paper Q we presented results bearing on the breaking of isospin using the same technique 
applied here with a selection of five different forms of potential. On examining the phase shifts produced 
outside the range of the fit it was noticed that the prediction was much better for two of the models, sums of 
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local Yukawa and exponential potentials for each partial wave. From fitting the data from 30 to 50 MeV these 
two models were able to predict the existence of the 33 resonance, and indeed the position of the resonance 
to within about 10% in kinetic energy. Thus it seemed reasonable to extend the analysis to somewhat higher 
energy with one (or both) of these two models. 

The Yukawa potential has the advantage that it naturally represents particle exchanges. As will be 
discussed in Sec. 3, the Klein-Gordon equation contains a term which is quadratic in the potential. For a 
Yukawa potential, such a term is singular (c>c l/?"^)- While solutions of this equation are obtainable, the 
result would seem to be more physical if a cutoff is introduced to soften the potential at small values of r 
(perhaps due to the intrinsic size the of quark-pion system). For an analysis employing particle exchange 
(hence approximately related to Yukawa potentials) see the recent work of Timnicrmans Q . 

One goal of this article is the representation of the pion-nucleon interaction in a simple and transparent 
manner. The exponential potential lends itself to that end without introducing a singular potential. Another 
advantage of the exponential form is that it might better represent the interaction of the pion with quark 
distributions within the nucleon. (The density calculated from a bound-state solution of a wave equation 
with a 1/r potential is exponential, as in the hydrogen atom and a three body system tends to follow the 
same density Q.) 

We have also made a short study utilizing the Gaussian potential, even though it is not considered to have 
a strong physical basis, to estimate which results are sensitive to the form of the potential (see Sec. VH). 

The article is organized as follows. Section 2 summarizes our choice of elastic-scattering data sets and 
briefly discusses their consistency. Section 3 reviews Jost's method for determining the s-wave amplitude 
for sums of exponential potentials. Our method of evaluation of the subthreshold amplitudes is described in 
Section 4. Section 5 begins by presenting numerical values of the potential parameters determined from a 
fit to the data. This is followed by subsections which present the consequences of the fit, such as the phase 
shifts, scattering lengths, coupling constant, sigma term, off-shell amplitudes, partial integrated elastic cross 
section, and polarization asymmetry. The results are summarized in Section 6. 



II. DATA 

The data sets that we have considered come from, for the most part, experiments dedicated to the mea- 
surement of pion-nucleon scattering. We have excluded experiments which have not been published or which 
had the measurement of pion-nucleon scattering as a secondary goal. The sets used were: 

Sigg The one atomic measurement by Sigg et al. is very important in determining the low-energy 
behavior of the s-wave amplitudes. It alone fixes, to large extent, the scattering length of the tt^ -p system. 
We found that predictions of this scattering length from our fits to scattering data were always near to their 
value and, when included, the fit adapted itself easily to the value and usually fit it with a very small x^- 
Thus we see no reason to question the validity of this point and have used it in all of the analyses discussed 
below. 

Brack The data from J. Brack et al. seem to be of high quality. The smallest errors are (in general) 
those quoted in this work, which contains 62 data points for scattering of 7r+ and tt~ at 29.4, 45.0, 66.8 and 
86.8 MeV. 
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Frank The data from J. Frank et al. g were the first modern data to be pubhshed contradicting the 
Bertin et al. |^ data. This work contains 162 data points for scattering of 7r+ and at 29.4, 49.5 and 69.6 
MeV. 

Auld The data of E. G. Auld et al. are somewhat old now and have large error bars but are generally 
consistent with the modern data. These data contain 11 points for scattering of 7r+ at 47.9 MeV. 

Ritchie The data from B. Ritchie et al. |jll| consist of 28 points of back-angle positive pion scattering at 
65.0, 72.5 and 80.0 MeV. We have found, like Timmermans that the normalization seems to be flawed in 
this data. We have floated the normalization of the lowest energy data set. (This data is from a secondary 
experiment.) 

Wiedner The data of U. Wiedner et al. |jl^ represent the flrst to appear in print from PSI. This set 
consists of 7r+ and 7r~ scattering at 54.3 MeV. Our fit indicates that the cross sections for the negative pion 
data are too high by about 14%. 

Joram I The first paper of Joram et al. presents data in the important Coulomb-nuclear interference 
region for 7r+ and 7r~ at 32.20 and 44.60 MeV. The set contains 80 data points and seems to be consistent 
with the rest of the data base. 

Joram II The data in the second paper of Joram et al. ||l^ present larger angle cross sections of 7r+ and 
IT" at 32.3, 44.6 and 68.6 MeV. The authors made single-energy fits to the data in the experimental paper. 
For the 7r+ at 32.3 MeV the best /N that they obtained, in combination with other data at this energy, 
was 121/58, while for the 7r+ at 44.6 MeV the best was 95.6/46. It was pointed out in Ref. ||l^ that the 
scattering lengths obtained from these single-energy fits lead to values of the ttNN coupling constant outside 
of an acceptable range. We have not included these two data sets in our fit. We also have dropped the 53.42 
degree point in the 7r+ data at 68.6 MeV which is completely out of line with the rest of the data. This 
leaves 32 data points in this set. 
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FIG. 1. Cross sections from tt proton scattering around 50 MeV. The Bertin data would agree with the prediction 
given by the dotted curve if they were consistent with the other data sets. The solid, long dash and dash-dot curves 
come from our fit and correspond to the Karlsruhe, Brack and Frank data sets. 



The data of Bertin et al. g| were not included in the general fits since they (or at least their normalizations) 
seem to be inconsistent with the modern data. In order to make contact with previous analyses and get 
some feeling for the impact of the data sets, we make fits using the Bertin data as the only 7r+ data and the 
full modern 7r~ data set described above. Figure 1 shows the tt"*" data around 50 MeV compared with one 
of our fits illustrating the apparent discrepancy. 



III. SOLUTIONS FOR THE AMPLITUDES: JOST REPRESENTATION 

At low energies and for each partial-wave, the tt"*" p elastic process is described by a single- channel Klein- 
Gordon (KG) equation. The tt" p elastic and charge-exchange scattering are described by a two-channel 
KG equation. That the effect of the (tt", 7) reaction on the hadronic channels may be ignored was justified 
in Ref. |^ where this procedure is discussed in detail. 

In fitting the data, we solve the KG equation by standard numerical procedures. The potential, V , is 
included in the KG equation through the substitution uj lo ~ V , where uj is the center of mass energy of 
the pion (actually the reduced energy in the fits); i.e. for both electrostatic and strong interactions, V is 
taken to be the fourth component of a Lorentz four-vector. The resulting equation is 
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{\/^ + - 2u;V + V^)iIj = 0. 



(1) 



where k is the center-of-mass momentum. 

However there are several calculations that wc would like to make which include only the strong interaction 
(the sigma term, the isovector s-wave amplitude at the CD point, and the ofF-shell amplitudes) for which 
the Coulomb force is not included. In this case it is very convenient to use (for the s-wave) the expressions 
developed by Jost p6[ |. 

We now give a brief outline of the method and its extension to the off-shell case. Our potential contains 
(at most) two terms of exponential form. Since the potential appears linearly and quadratically in the KG 
equation, a two-term basic potential leads to an effective Schrodinger equation with a 5-term potential. 

Consider the solution of the Schrodinger equation for a sum of exponential potentials: 

N 

F(r)=^A,e-'^^'-. (2) 

Jost writes the solution for the s-wave, f{k,r), as 

/(fc,r)=e-*^-'-^C„(fc)e-™"'- (3) 

a 

where the subscript a is a compound quantity consisting of a set of N integers. For example, for a three-term 
potential 

a=[j,k,l], J, k, I = 0,1, 2... (4) 
The coefficients Ca{k) are given by the recursion relation 

"''"''^^ ' ma{ma + 2ik) ^ ' 

where 

nia = myj,kj] = jMi + ki-i2 + llJ-z- (6) 

The recursion is started with 

C[o,o,o] = 1, Cy^i.kA = (^[j-iA = C\j^k-i\ = (7) 

and is built up by first computing all coefficients with the sum of indices equal to one, then two, etc. with 
no negative index. 

The solution with the proper boundary conditions at the origin with an incoming spherical wave and with 
unit amplitude at infinity is 

where 
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/(-fc,o)- 



(9) 



These expressions can be used to calculate the values of the S-matrix for any value of k. We shall be 
interested in purely imaginary values for the calculation of the s-wave contribution to the E-term. 

In order to calculate the off-shell amplitude we consider the wave function for real (positive) values of k. 
In this case we can write (for real ) 



f{-k,r) = f*{k,r) with S = e^'^ and f{k,0)=e'^(3 
where (3 is real and positive. We can now write the wave function (Eq. H) as 



■)p{k,r) = - 



2ikr 



The s-wave off-shell amplitude is defined by 

Fo{k,q) = j drr'^jo{qr)V{r)i){k,r) 

and wc find 



Fo{k,q) = -fT^'^ 



A. 



{nia + Mj)^ + g2 - fc2 -I- 2ik{ma + fJ,j) 



This function is evaluated in Sec. VE 



(10) 



(11) 



(12) 



(13) 



IV. SUBTHRESHOLD EXTRAPOLATION 



The subthreshold regime, the region below the elastic threshold, is of interest in studies of chiral-symmetry 
breaking, a measure of which is the S term. On-shell subthreshold amplitudes have been evaluated over the 
years by the use of dispersion relations It is interesting to compare those results with the values 

given by the KG equation. 



A. Dispersion Relations 

Because the "experimental" S-term is defined at the unphysical CD point (s = m^, t = 2/i^), any 
determination of it is, to a degree, model or theory dependent, s-channel and w-channel data exists only for 
negative t, therefore fixed-i dispersion relations must either rely on the rapid convergence of the partial-wave 
series outside of the physical region at t = 2^^ as in Cheng and Dashen's method or must use additional 
techniques (such as dispersion relations at fixed v) to extend the amplitudes to t — 2/i^. This extrapolation 
is complicated by the presence of the crossed reaction tttt —^ NN which produces a nearby branch cut in 
the t-channel beginning at t — . Some knowledge of the effect of this cut is obtained indirectly from vrvr 
elastic scattering, but this analysis is, to an extent, model dependent. Some formulations of the dispersion 
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relation approach rely heavily on very accurate knowledge of the p-wave scattering volume [Ell when, in 



fact, the S term is dominated by the s-wave amplitude [p2 23 



Interior Dispersion Relations (IDR) |2C|] and Hyperbolic Dispersion Relations have also been used to 
determine the E-term. In these cases, the curves along which the dispersion relations are written may be 
chosen to pass directly through the CD point. However, these dispersion relations involve integrals over 
the entire t -channel cut, a long portion of which is unphysical < ^ < 4m^). In this case, the entire 
i-channel dispersion integral (the "Discrepancy Function" ) must be determined and extrapolated to the CD 
point. Thus, even in the classic dispersion-theory determination of the S-term, some model dependence is 
present. 



80 




u (MeV) 

FIG. 2. Contributions to the E term for the s-wave (dash-dot curve) the s- and p-waves (dashed curve) and the 
fuU amplitude (solid curve) from the polynomial expansion from Ref. [M From this graph the p-wave is estimated 
to give a correction of ~ +1 MeV. The sum of all of the non-s-wave contributions gives ~ -1.5 MeV. 

It is instructive to compare the relative contributions of the s and p waves to the non-flip pole-subtracted 
amplitude which is real in this region and, as we shall show in the next section, is proportional to E 
at the CD point. To this end we have performed the partial-wave projections of as computed from 
the subthreshold polynomial expansion in v = {s — u)/Am and t from Hohler et al. j2^. For purposes of 
illustration we plot vs. v along a curve of constant path parameter a = — 0(s,t)/t^, where 4>(s^i) = 
t \su — (vr? — /i^)^] is the Kibble function, which vanishes on the boundary of the physical region. The 
parameter, a, is fixed at the value — -|- /^^/2. This path has the advantage that it passes through both 
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the s-channel threshold {ly = fi) and the CD point (at which 1^ = 0). It also arrives at the CD point tangent 
to a curve of constant t as was specified in the paper of Cheng and Dashen. (This is the path used in IDR 
calculations.) For comparison, the s and s+p wave partial-wave contributions to are plotted along with 
the full amplitude in Fig. ^. The amplitudes have been multiplied by a factor so that they reduce to the 
S-term at = 0. It is seen that the s wave contribution is similar to the full amplitude between the CD 
point and threshold, becoming quite close at the CD point. Notice that the s and s+p contributions cross 
slightly below the CD point, where cos 6* vanishes. This zero in cos 6* suppresses the p-wave contribution at 
the CD point p2] . The value of cos 6 at the CD point is 



1 

cos6l = f . (14) 



The value of t for which cos 9 is zero (for v — 0) is approximately 



which differs from the value at the CD point by about 1/2 of one percent. 



B. Klein-Gordon Approach 

In the present approach the subthreshold dependence is derived directly from the pion-nucleon data itself. 
The question of the validity of the function extracted then depends on the accuracy of the data and the 
suitability of the form used for the potential. 

While the amplitudes produced from the KG model do not automatically contain the full analyticity and 
crossing properties of the invariant amplitudes, much of the same physics can be included by considering 
that the source of the potential is the i-channel cut (tttt scattering and tttt — > NN). The potential can be 
written in terms of an integral over the discontinuity along the t axis as done in the case (for example) of 
the NN potential [|5) 

V{r) = / Pit) dt. (16) 

Since a Yukawa potential can be written as an integral over an exponential form 



oc 

xr 



e-^'^'dx (17) 



the ttN potential can be expressed as a sum of exponential potentials. 
From the definition of the non-spin-flip amplitude, 

G{s,t) = ^[(2f + 1)/,+ +£/,_]P,(a:), (18) 

and the projection of these partial waves from the invariant amplitudes A{s, t) and i?(s, t) (see Ref. Eq. 
2.7 for example) we can write 

167rsG(+' (s, t) = [{W + m)2 - ^2] [^(+) (^^ t) -\- [W - m)B^+^ [s, t)] 
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+ x[iW ~ mf - fi^][-A^+\s,t) + {W + m)B'^+'>{s,t)] (19) 

where W = ^/s and x {— cos 6) is a function of s and t. The same result is available from Eq. A 2.32 of 
Ref. [|l^. We omit the superscript "(+)" for the next few equations for clarity. For t ~ 2/i^ the Born term 
contribution to B, ( — i ^ — |, becomes 
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If we define B{s, t) as the _B-amplitude with the pole term removed then 



167rsG(s, 2/x2) = [{W + mf ~ ^?] 



A{s, 2i?) + {W - m)i3{s, 2n^) 



W + m 



+ x[{W - mf - /i^l 



J2l_ 

W — m 



(20) 



A{s, 2^?) + {W - m)B{s, 2fi^) 



29' 



W + m 



16nsG{s,2fi^ 



2g^fi^x 
W — m 



x{W - mf[~A{s, 2i?) + (VF + to)B(s, 2^?)] + x^?[~A{s, 2^^) + {W + m)B{s, 2^'^)] 



(21) 



Now define G(s, 2/i^) = G{s, 2/i^) — STTs{w-m) remove the pole from the Pn partial wave. The subtracted 
pole term should be evaluated with the fitted position and residue so that the pole is exactly removed from 
the Pii partial wave of G{s, 2/i^). Then, at the CD point (s — m?), 



{An? - ^JL^) 



o2n 



A{m^,2^?) 



= l&TTm^G{m^, 2fi^) + xcDfJ-'^ -A{m^,2i?) + 2mB{m^, 2^?) 



(22) 



The last term is negligible at the CD point as we will show shortly. Thus, neglecting the factor (1 — j^); 



A^T&+\m\2^?) = A^+\m\2^?) - ^- = E//^ 

m 



(23) 



where /^r = 93.2 MeV is the pion decay constant. The equation corresponding to the second equal sign has 
been given many places (see for example). The corrections to this expression of the order (^)^ psf . 
The ratio of the last two terms in Eq. ^ is approximately 

xcD^Ji^[-A + 2mS]/^/(WS) « 0.002 

where we have used A « g'^/m w 191 GeV"\ 2m]3 « -311 GcV"\ E « 0.060 GeV, and xcd ~ - fi'^ / Am'^ . 
Thus we may safely ignore the last term in Eq. (22). 

Thus, to a good approximation, the sigma term can then be evaluated as S = 47r/^G'^+-' where 
is evaluated at the CD point. We approximate G'+' in terms of s-channel isospin amplitudes: G^+^ = 
i(GU2Gi). 

2 

The CD point occurs at center-of-mass total pion energy, a; = ^ « 0; at this point, the effective potential 
U = 2loV — V"^ ~ —V^ is attractive whatever the sign of V . Above threshold the cross term dominates, 
and U may be either positive or negative. We are particularly interested in the s-wave scattering amplitudes 
because, as discussed in the previous section, at the CD point cos « and the p-wave contributions to the 
sigma term are suppressed. At threshold one of the s-wave amplitudes, 6*1/2 is positive, and the other 63/2 
is negative. We define renormalized amplitudes by 



= 47r/2G^+' « ATiflSi « 553 Si 



(24) 
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where Sj are the s-wave scattering ampHtudes in fm. and / = 5, §• When evaluated at the CD point 
S = ^(Sj. + 2E3) is the sigma term. 

As u! passes below threshold (w = /i) these two amplitudes will keep their signs (they are real in this 
region) but ai uj = they must be of the same sign (and positive) because the square term now dominates 
(unless the squared potential is strong enough that there is a bound state). 

In the low-energy physical region there is a correlation between the range and strength of each potential; 
i.e. two fits are, to lowest order, comparable if the volume integrals of the potentials are equal. However, in 
the subthreshold region near the CD point the V'^ term dominates 2ujV (w « 0); the correlation no longer 
holds. Therefore to determine the amplitude near the CD point with our model, the strength and range of 
the potentials must be independently determined. We believe that the available data is sufficiently accurate 
that this determination can be made, at least approximately. 

Since A^~^(m^, 2/i^) = and the Born term is zero at the CD point, we find 

G(-'(™^ V) ~ 0, (25) 

where the zero is of order (j^)'^ times a typical amplitude. Since the isovector amplitude G^^\m'^ ,2^'^) 
vanishes, it follows that the isovector combination — Sj. — E3 should vanish (assuming s-wave dominance) 
at the CD point so that the two amplitudes defined in Eq. ^ should cross at that point. This vanishing 
is not automatic in the KG model, which doesn't have s <-> u crossing built in. In performing a number of 
fits without this requirement we found that Ei and E3 do indeed naturally cross near v = Q. However, the 
isovector amplitude varied between 10 and 100 MeV at the CD point. Even so, we obtained stable values 
for the scattering lengths while values of the E term varied from 33 to 100 MeV. 



Using the Jost technique described in Section III we evaluated S at the CD point at each step of the fit. 
By including a contribution of 

/ E 3 - E , \ 2 

(26) 



to the it was possible to force the vanishing of I]^(i^ = 0) to a reasonable degree of accuracy. The 
result is that the value of the S term obtained is much more stable. See Ref. for a related discussion of 
subthreshold constraints on potential models. 



Partial Wave 


ai (MeV/c) 


Ai (MeV) 


02 (MeV/c) 


A2 (MeV) 


Ri (fm) 


R2 (fm) 


S3 


723.3 


901.0 


485.3 


-41.4 


0.945 


1.408 


Si 


750.8 


-1.6 


310.5 


-58.0 


0.910 


2.201 


P33 


891.7 


399.2 


486.7 


-1548.6 


0.767 


1.404 


Pl3 


813.9 


-5584.2 


528.5 


834.3 


0.840 


1.293 


P31 


631.4 


1800.5 






1.083 




Pll 


720.3 


-3540.3 






0.949 





TABLE I. Ranges and strengths for the potential obtained from a fit to the data. The radii are calculated in the 
rms sense. 
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FIG. 6. P33 wave phase shifts The curves and symbols have the same meaning as in Fig. 



V. RESULTS OF THE FIT 



The form of the potential for each partial wave was taken as 



(27) 



The second term was not needed for two of the partial waves. As an example of the values obtained, the 
strengths and ranges are given in Table | for the case 11 in Table 2. We see that the ranges corresponding to 
the dominant strengths mostly have values from 500-900 MeV/c. An exception is the partial wave Su where 
the value is ~ 310 MeV/c, barely above 2 pion masses, the minimum acceptable value without creating an 
anomalous threshold. That such a long (spatial) range is needed to fit this partial wave was noted by Ref. 



We have tested individual data sets to see their influence on the fits. The results presented in Tables O 



and II] show the effect of adding, one by one, the data sets in the order shown. Below the solid lines are fits 



to single groups of data. The second line that appears in most entries of Table [II corresponds to a second 
fit which gives an idea of the uniqueness of the fits. 



Brack 



X /n (data) 
64.66/62 



X /n (norm) 
9.88/10 



S (MeV) 
42.86 



a-i (/i ^) 
-0.083 



ai 



(A*-') 
0.174 



f 
0.0759 



16 



+Frank 


231.78/228 


21.85/16 


48.43 


-0.085 


0.174 


0.0762 


+Auld 


243.28/239 


23.20/17 


48.59 


-0.085 


0.174 


0.0762 


+Ritchie 


274.58/267 


27.68/20 


48.66 


-0.085 


0.174 


0.0762 


+Wiedner 


316.27/306 


33.09/22 


50.01 


-0.085 


0.174 


0.0762 


+Joram I 


432.67/386 


35.42/26 


48.90 


-0.081 


0.172 


0.0751 


+Joram II 


497.93/417 


42.01/30 


53.59 


-0.082 


0.172 


0.0753 


Frank 


144.15/166 


4.56/6 


48.34 


-0.083 


0.173 


0.0757 


PSI 


153.49/150 


8.33/10 


67.45 


-0.095 


0.180 


0.0793 


Bertin+TT^ 


278.56/243 


15.96/18 


52.92 


-0.105 


0.183 


0.0817 



TABLE II. Results of the study, is separated into contributions from the individual data points and from the 
experimental normalization uncertainty. The fits were made including a contribution to the total from Eq. ^ of 
AiSs = 10 MeV. The line labelled PSI includes both the Wiedner and Joram data. 







xVN (Data) 


xVN (Norm) 


E (MeV) 


as (m ^) 


ai (/i ^) 


f 


1 


Brack 


66.37/62 


9.70/10 


50.51 


-0.082 


0.174 


0.0757 


2 




65.09/62 


10.52/10 


45.53 


-0.083 


0.174 


0.0759 


3 


+Frank 


234.69/228 


21.56/16 


50.77 


-0.085 


0.175 


0.0764 


4 




230.63/228 


22.63/16 


46.32 


-0.086 


0.174 


0.0764 


5 


+Auld 


245.35/239 


22.88/17 


50.76 


-0.085 


0.175 


0.0770 


6 




242.43/239 


23.46/17 


46.32 


-0.086 


0.175 


0.0766 


7 


+ Ritchie 


276.50/267 


28.21/20 


50.79 


-0.086 


0.175 


0.0766 


8 




273.59/267 


27.94/20 


45.36 


-0.086 


0.174 


0.0764 


9 


+Wiedner 


317.40/306 


35.12/22 


49.86 


-0.086 


0.176 


0.0768 


10 




315.21/306 


33.12/22 


45.53 


-0.087 


0.175 


0.0768 


11 


+Joram I 


431.89/386 


37.97/26 


48.98 


-0.083 


0.174 


0.0759 


12 




432.83/386 


35.58/26 


45.14 


-0.083 


0.173 


0.0757 


13 


+Joram II 


495.96/417 


44.57/30 


49.22 


-0.084 


0.174 


0.0761 


14 




500.67/417 


42.64/30 


45.99 


-0.085 


0.173 


0.0761 


15 


Frank 


139.07/166 


4.41/6 


48.67 


-0.083 


0.174 


0.0757 


16 




139.27/166 


4.71/6 


46.51 


-0.083 


0.174 


0.0757 


17 


PSI 


153.32/150 


8.38/10 


65.17 


-0.095 


0.180 


0.0793 


18 


Bertin-f-TT^ 


271.40/243 


17.47/18 


50.68 


-0.103 


0.184 


0.0826 


19 




273.53/243 


15.87/18 


48.94 


-0.105 


0.184 


0.0820 



TABLE III. Results of the study. The fits were made including a contribution to the total x^ from Eq. ^ of 
A_Be = 1. 
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FIG. 7. S-wave amplitudes for pion-nucleon charge exchange, measured (solid lines) and predicted by isospin from 
the present fit (dotted lines) 



A. Phase Shifts and Low-energy Parameters 



Figures ||, U, || and ^ show the phase shifts resulting from the fit 14 of Table VU. 
We obtain for the isoscalar and isovector scattering lengths: 

6o = 0.001 ±0.003 and 6i = -0.086 ±0.003 (28) 

in reasonable agreement with the values obtained by Sigg et al. |5j. 

A recent QCD calculation by Bernard et al. |3lf| gives -0.096//"^ <bi< -0.088Ai"i with which we are m 
marginal agreement. 

The Chew-Low theory |3^ predicts the values of the scattering volumes for all of the p-waves. For the 
P33 the prediction is 0.191/j,~'^ (using = 0.0765), in reasonable agreement with the values obtained. 

An interesting feature of the P31 and P13 partial wave phase shifts (Figure U) is that they are very nearly 
equal at low energies, becoming essentially identical at threshold. Above 20 MeV in kinetic energy the curves 
separate. The Chew-Low theory predicts that these partial waves are always equal. Improved models [ |3^ , |3^ 
give different magnitudes for these scattering volumes, but the symmetry of the hamiltonian requires that 
they be equal. 

It is not clear how fundamental this symmetry should be considered to be, and thus on what level we should 
expect it to hold. It has been known from some time that this spin-isospin symmetry holds in the Skyrmion 
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model . Recently it has been shown |Q that the relations among partial waves of the Skyrmion model 
hold more generally in the large l/N^ expansion of QCD. Since the predictions are for infinitely massive 
nucleons the symmetry might well be expected to be broken above threshold where recoil effects become 
important. 

Comparisons with data have always shown this symmetry not to be satisfied in nature, but we find that 
it is satisfied within the errors. We find for both scattering volumes —0.028 ± 0.003/1^''. 

Figure H shows that the KH80 phase shifts follow our curve for values above 50 MeV but deviate at 
20, 30, 40 and 50 MeV, just at the Bertin energies. Table VA shows that, indeed, the analysis using the 
Bertin data results in a P31 scattering volume considerably larger than the P13. 

For the P33 phase shift (Fig. |^) we note that the data fit in the region of 30 MeV to 86 MeV give the 
behavior of this resonant phase shift with a reasonable accuracy. 

The Pii phase shift (Fig. ^ is seen to have a negative excursion at low energies and to cross zero around 
140 MeV, as has been typical in previous fits. However, in the present case, it is the low-energy data alone 
which find this behavior, rather remarkable in view of the small size of the amplitude. 







ba 


bi 


P33 


Pl3 


P31 


Pii 


1 


Brack 


-0.050 


-0.068 


0.181 


-0.031 


-0.029 


-0.102 


2 




-0.049 


-0.064 


0.178 


-0.031 


-0.031 


-0.098 


3 


-l-Frank 


-0.050 


-0.068 


0.179 


-0.027 


-0.028 


-0.104 


4 




-0.049 


-0.065 


0.176 


-0.027 


-0.029 


-0.102 


5 


+Auld 


-0.050 


-0.068 


0.179 


-0.027 


-0.027 


-0.105 


6 




-0.049 


-0.065 


0.176 


-0.027 


-0.028 


-0.102 


7 


-l-Ritchie 


-0.050 


-0.068 


0.180 


-0.027 


-0.026 


-0.106 


8 




-0.048 


-0.064 


0.177 


-0.028 


-0.028 


-0.102 


9 


+Wiedner 


-0.050 


-0.073 


0.178 


-0.026 


-0.029 


-0.108 


10 




-0.048 


-0.064 


0.177 


-0.027 


-0.028 


-0.102 


11 


+Joram I 


-0.049 


-0.071 


0.179 


-0.027 


-0.028 


-0.108 


12 




-0.048 


-0.063 


0.177 


-0.028 


-0.028 


-0.103 


13 


+Joram II 


-0.049 


-0.071 


0.178 


-0.025 


-0.028 


-0.106 


14 




-0.049 


-0.064 


0.176 


-0.026 


-0.027 


-0.102 


15 


Frank 


-0.049 


-0.071 


0.169 


-0.015 


-0.025 


-0.109 


16 




-0.049 


-0.065 


0.168 


-0.014 


-0.025 


-0.108 


17 


PSI 


-0.055 


-0.083 


0.183 


-0.027 


-fO.392 


-0.119 


18 


Bertin+7r~ 


-0.049 


-0.076 


0.178 


-0.023 


-0.033 


-0.100 


19 




-0.049 


-0.070 


0.177 


-0.024 


-0.034 


-0.094 



TABLE IV. Effective ranges (in units /i ^) for the s waves (ba and bi) and scattering volumes (in units /i from 
the fit. The entries in this table correspond to those in Table pTD. 
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B. Isospin Breaking 



We have calculated the prediction for the charge- exchange amplitude for each of the fits shown in Table 



III in the same manner as discussed in Ref. 



The charge-exchange data fit are the same as before |38|-|40(|. Note that the fit to this data is direct, 
it requires no model and could have been done with a simple polynomial. The normalization of the data 



is constrained with the use of the TRIUMF data |39 4Q]. While these data provide only limited angular 



information, they give a relatively accurate value for the integrated charge-exchange cross section and thus 
provide an important normalization constraint. 

While in the previous work we demonstrated a small dependence on the model used for the charged 
pion scattering, here we show (Fig. ^ a comparison of the charge-exchange amplitudes determined directly 
from the measurements (the same as in with the predictions of six of the fits from Table III (1, 2, f 1, 
12, 15 and 16). As can be seen, the variations among the data sets is small compared with the discrepancy. 

It is worthwhile to discuss the size of the result. It is possible to obtain a large fraction for the breaking 
because of the smallness of the basic amplitude itself. Let us compare with the NN case. If we assume that 
the potential causing the breaking is the same for ttN and NN (as happens to be the case for pcj mixing |^ ) 
then (in Born approximation) the breaking amplitude should be smaller for the pion-nucleon case by a factor 
of the ratio of the pion mass to the nucleon mass (about 0.14). However, a typical pion-nucleon amplitude 
at low energy (0.1-0.2 fm) is a factor of 100 smaller than a typical nucleon-nucleon amplitude (15-20 fm). 
Hence the percentage of the breaking should be around 10 times larger in the pion case compared to the 
nucleon-nucleon case. The breaking amplitude observed here ('^ 0.012 fm) would correspond to a breaking 
in NN case of around 0.1 fm which is too small to be observed at the present time. 

The predictions with the Bertin data, from line 19 in Table HI (not shown in the figure) agree with the 
measured charge exchange (as previously observed [|l[). 

A possible explanation has been advanced by Piekarewicz ]4^ ] for this breaking in terms of the difference 
of vr" pion-nucleon coupling constants for the neutron and proton. This explanation would put the breaking 
entirely in the charge- exchange channel or the ay coefficient in the notation of Ref. 1 43 1 . 



C. The Pion-nucleon Coupling Constant 

The GMO (Goldberger-Miyazawa-Oehme) jJij sum rule provides a useful relation between the isovec- 
tor combination of the ttN scattering lengths and the ttNN coupling constant This sum rule is obtained 
from a forward dispersion relation for the invariant amplitude subtracted at lolab = and evaluated at 

UJLAB = M ill- 



4.^(ar-a3) = ^^+4.J, (29) 
3m/i — Ljg 



where 

1 r a-{k) 



27r2 7o ijj(k) 



dk. (30) 
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In this expression k is the incident pion laboratory momentum and ujb = — /i^/2m. The isospin-odd total 
cross section is defined by a~ = 5(0'- — c+) where a± = axiT^^p)- 

To evaluate J we have used the SM95 |46| phase shift analysis for fc < 2 GeV/c, our own parameterization 
for 2 GeV/c < k < 4 GeV/c (a fit to total cross section data taken from the Review of Particle Properties 
p7t), and a Regge fit for 4 GeV/c < k < 240 GeV/c taken from the same source. The integral was truncated 
above 240 GeV/c. The resulting value of the integral is J=(-1.308+0.068+0.157) mb =-1.081 ± 0.005 mb 
where the contributions from each momentum interval is shown. If we were to assume that the Regge fit is 
valid to infinite momenta, the contribution to J coming from above 240 GeV/c is 0.030 mb, which would 
lead to J=-1.051 mb. For the following discussion we will use the truncated value, J=-1.081 mb. This value 
can be compared with J=-1.072 mb obtained by Ref. |Q and J=-1.077 mb quoted in Ref. |Q from an 
unpublished preprint by R. Koch. Thus it would appear that uncertainty in the determination of J is the 
order of 1/2 to 1%. Since it contributes about 1/3 of the value of the error in from this source is 
less than 1/3 %. Locher and Sainio Q concluded, however, that the uncertainty in J was slightly larger, 
leading to a 1% error in 

Thus, we have the simple formula for (ai and 03 in pion mass units) 

/2 = 0.0269 + 0.1904(ai - as). (31) 
Some of the results of our fits to the low-energy n^p elastic scattering data are shown in Tables |l^ and 



inl. The resulting average value of the coupling constant is = 0.0764 ± 0.0007 where the error quoted 
includes our fitting error only. 

The values for the scattering lengths advocated by SM95, 03 = — 0.087/i~^ and ai = 0.175 lead 
to — 0.0768 when used in our relation. Arndt et al. ^ quote a value of p = 0.076 ± 0.001 while 
Markopoulou-Kalamara and Bugg (50) found = 0.0771 ± 0.0014. Timmermans quotes a preliminary value 
from his ttN analysis § of 0.0741 ± 0.0008 (statistical error only). 

Thus recent analyses of the charged pion coupling constant from pion-nucleon scattering seem to be in 
moderately good agreement. It is very interesting to know if this value is consistent with that obtained 
from the nucleon-nucleon interaction since that comparison serves as a check on our theory of the strong 
interaction. 

The Nijmegen group [|l| found p = 0.0748 ± 0.0003 for the charged pion coupling constant and = 
0.0745 ± 0.0006 for the neutral pion. Ericson et al. |5^ find from the analysis of np charge-exchange data, 
/2 = 0.0808 ±0.0017. 

While results here would seem to confirm (perhaps even with smaller errors due to the determination 
of the scattering lengths directly from the low-energy data) the results of other analyses of pion-nucleon 
scattering, there is an important caveat which should be mentioned in regard to all of the analyses with 
the GMO sum rule. It is the isovector scattering length which enters in the GMO relation. As we saw in 
the previous section, there is a strong indication of an isospin breaking from the comparison of the elastic 
scattering determination of this amplitude with its determination from charge exchange. If the explanation 
advanced in Ref. |Q is correct (so that all of the breaking is in the charge- exchange channel) then the 
GMO relation will lead to a correct value of If there is some breaking in the elastic scattering channels 
then that correction should be made before applying the GMO relation. In the extreme case that all of the 
breaking comes in the elastic amplitudes (an amplitude of the form of p3[ , possibly from p — uj mixing. 
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with the opposite sign than predicted) and the isovector amphtude from charge exchange is the correct one 
then the pion-nucleon coupUng constant would return to the "large" values. It is not possible to reliably 
determine (from GMO) until the question of the origin of the isospin breaking is resolved. 



D. The S Term 



Figure g shows the amplitude defined by Eq. |2J as a function of v (which is roughly equal to the center of 
mass pion energy, w, in this region). Since above threshold the amplitude, the real part was used. The result 
of the solid curve at v = represents the s-wave contribution to S in our model. The solid and dash-dot 



curves were calculated using the Jost solutions for fit 14 in Table [II . The kinematic singularity at threshold 
is clearly visible. 
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FIG. 8. S-wave E term amplitude vs f from the present work (solid curve) compared with the pole-subtracted 
expansions of Nielsen and Oades ^ (short dash) and Hohler ||| (long dash). The two isospin components from 
our work are shown with the dash-dot curves. 



While this plot is given as a function of u, the more relevant variable for extrapolation may be the center 
of mass momentum. Since the CD point lies at fc « i/x, a circle in the complex plane with radius passes 
through the real axis at a kinetic energy around 57 MeV, or in the center of range studied in this work. 

The values obtained for E lie around 48±4 MeV. The error quoted is determined from the variation among 
data sets and does not include the model error. It is interesting to note that there is very little difference 
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between the results from the new 7r+ data and that of Bcrtin ct al. Q . This is perhaps understandable from 
the observation that the Bertin data for the cross section are 20% higher than the prediction from the fits 
(see Fig. 1) so the 7r+ p amplitude can be expected to be 10% higher. Since the tt"*" and tt^ contribute 
equally to S, the difference of the analysis between the two data groups can be expected to be of the order 
of 5% or 2-3 MeV. Indeed the difference observed is of that order. 

While the usually accepted value of S is around 65 MeV this is not the first time that a smaller value 
has been obtained. Ericson found a similar value (44 ±6 MeV) |^^. The modern approach seems to 

be to consider this value as not including the part of the S term coming from the cut in t. This additional 
part (about 11 MeV) is added on after the extrapolation of the pion-nucleon amplitude. If we compare our 
value of S with of Ref. |^l| (before correction for the tttt channel) we are in good agreement. However, 
the TTTT channel is implicitly included in the present analysis since an integral over the discontinuity of the t 
cut gives the potential. Thus the singularities in t are present in the model and we must consider the value 
obtained as including this correction. 



E. Off-shell Amplitudes 
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FIG. 9. S-wave off-shell amplitudes. The solid line corresponds to the isoscalar combination of the two amplitudes. 



23 



1.1 
1.0 
0.9 
0.8 

0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 
=0.1 









- 










\\ 

\ 








- \ 

\\ 

\\ 


^_^P.3^^^ 












\ \ 
p < 

31 \ 


P33 















400 



800 

q {MeV, 



1200 
'c) 



1600 



2000 



FIG. 10. P-wave ofT-shell amplitudes. Plotted is the ratio of the off-shell value to the on-shell value at 5 MeV. 
A knowledge of the off-shell amplitude is necessary for the treatment of pion-nucleus scattering. An early 



determination of these amplitudes was made by Landau and Tabakin 1 53 1 using the assumption of a separable 
potential and deducing the form of the off-shell amplitude from the energy dependence of the on-shell 
amplitude. Since the knowledge of the scattering amplitude was needed in the region where the scattering 
becomes inelastic, the potentials obtained were complex, not a very satisfactory situation. Londergan, McVoy 
and Moniz |5j] were able to find real functions for the off-shell dependence by considering only the potentials 
in the elastic channel, including the inelasticity by coupling to one other 2-body channel. 

In our case, since we are fitting to data only in the elastic region, there is no problem of a complex 



extension. We can calculate the off-shell amplitudes directly by the Jost formula given in Section |II| for 
the s waves or by direct numerical solution for the p waves. Figures ^ and ^ show the ratio of the off-shell 
amplitudes to the value on-shell at 5 MeV. 

Of particular interest is the isoscalar combination of the s-wave amplitudes. At low energy this combination 
nearly vanishes on shell, but this cancellation does not necessarily occur off shell. As can be seen from Figure 
^, for larger values of momentum there is no cancelation at all. 

Recent measurements have resulted in accurate threshold cross sections for vr" production in pp 
collisions. The most common calculation of pion production achieves the needed momentum sharing among 
the nucleons by rescattering of the meson after it is emitted from one of the protons. Since the s-wave 
7r°N scattering length is near zero, this contribution would seem to be small at threshold. Based on this 
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assumption, Horowitz et al. and Lee and Riska constructed models based on heavy-meson exchange 
which were able to explain the data. 

Using the fact that the tt" is produced "ofF-shell," Hernandez and Oset were able to explain the cross 
section using an estimated dependence of the off-shell behavior. Their results are uncertain, however, due 
to a lack of knowledge of this dependence. 

For the p waves (Fig. 10) it can be seen that the off-shell dependence of the P31 and P13 (as well as the 
P33) amplitudes is nearly identical below 300 MeV/c. This fact indicates that the long-range part of the 
interaction is the same for these three waves. 



F. Partial Total Cross Sections 



Recently transmission measurements have been made [q9|,|0| to determine the integrated elastic cross 
section beyond a fixed angle: "partial total cross sections" . While many of these measurements were made 
at higher energies than those treated here, there are several points in the energy region of interest. We 
have not included these points in the fit but now compare with the values obtained. Until recently the two 
measurements have differed but now there seems to be general agreement between them, but disagreement 
with the Brack data above 50 MeV. Our values are shown in Table 0. 



(MeV) 


Present Work 


Ref [|| 


Rcf @ 


39.8 


7.6 




8.5±0.7 


44.7 


8.9 




9.2±0.8 


45.0 


9.0 


10.8± 1.0 




51.7 


11.3 




11.8±0.8 


52.1 


11.5 


12.4± 1.0 




54.8 


12.6 




13. 2± 0.5 


59.3 


14.7 




15.8± 0.4 


63.1 


16.6 


18.0± 0.6 




66.3 


18.4 




20.4± 0.4 


67.5 


19.2 


20.7± 0.6 




71.5 


21.7 


23.8± 0.6 




80.0 


28.1 




29.7± 0.7 


92.5 


39.6 


43.3± 1.5 





TABLE V. Partial total cross sections in mb at the 30° limit. 
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G. Polarization Asymmetry 
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FIG. 11. Comparison of the prediction of the fit with the polarization asymmetry data of Wieser et al. [ ]6l|j6^ 

While there is no published data on the polarization asymmetry in low-energy charged pion scattering 
there is one measurement at 68.34 MeV which has appeared in conference proceedings |pl| and in a thesis 




32[. The prediction of our fits are in excellent agreement with these data as shown in Fig. 11 



H. Gaussian Fits 

While a Gaussian potential has little theoretical justification, some fits were made with this form to test 
the sensitivity to the model potential. If a result is stable under this change it may well be believed to be 



weakly dependent of the form of the potential function. The results are shown in Table VI. It is seen that 
the values for the pion-nucleon coupling constant remain within a very narrow range, the same range as for 
the exponential potential, and thus are not sensitive to the potential form chosen. 

xVN (Data) x'/N (N) S (MeV) ag (m"') (m"') f P^^^ P^^^ 

1 Brack 68.25/62 13.45/10 22.59 -0.085 0.174 0.0762 0.168 -0.029 -0.028 -0.061 

2 +Frank 231.72/228 26.64/16 22.86 -0.089 0.175 0.0771 0.166 -0.026 -0.026 -0.063 

3 +Auld 242.71/239 27.48/17 22.81 -0.089 0.175 0.0771 0.166 -0.026 -0.026 -0.063 
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4 +Ritchie 272.89/267 31.90/20 22.83 -0.089 0.175 0.0771 0.167 -0.026 -0.025 -0.063 

5 +Wiedner 314.50/306 37.36/22 22.89 -0.089 0.175 0.0771 0.166 -0.026 -0.025 -0.063 

6 -l-Joram I 441.25/386 39.86/26 22.86 -0.086 0.173 0.0762 0.167 -0.027 -0.025 -0.064 

7 +.I()ram II ."•)l().:iO/ 117 Ifi.OS/.'iO 23.11 -0.087 0.17:! 0.070 1 O.lfifi -0.02:) -0.025 -O.OCil 

TABLE VI. Results of the study with a Gaussian potential with the condition that A£e = 1 MeV. 



Another result which is only weakly dependent of the potential form are the scattering volumes of the P31 
and P13 partial waves. They are again found to be equal, albeit slightly smaller than for the exponential 
case. 

The values of other scattering volumes are altered significantly, especially the Pn wave which is small in 
this region and hence difficult to determine. 

The value of the S term is found to vary by less than ±1/2 MeV over the data sets, but with a different 
value than for the exponential potential. Thus, in the KG model, the extrapolation to the CD point depends 
on the form of the assumed potential. 
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VI. SUMMARY 



We have analysed recent low-energy pion-proton elastic scattering data. The principal findings are: 

1. The large bulk of the modern elastic scattering data is internally consistent and inconsistent with the 
older "Bertin" data. 

2. The observation of the isospin violation previously indicated was confirmed for the variations in the 
data set fit over a larger range in energy. 

3. The value of the subthreshold parameter, the S term, has been extracted with the exponential potential 
and the result is a value around 50 MeV, lower than previous estimates. The smallness of this value 
was shown not to be because of the change in data (which had little effect) but due to the model 
which was used. This conclusion is reinforced by the study with the Gaussian potential (considered 
unrealistic) in which a yet smaller value was found. 

4. The pion-nucleon coupling constant was extracted from the scattering lengths with the use of the 
GMO sum rule and a value (/^ = 0.0764 ± 0.0007) in agreement with most (but not all) modern 
determinations was obtained. 

5. Scattering volumes were extracted for the p waves. The parameters for the P13 and P31 were found 
to be the same within errors, as predicted from the Chew-Low theory and various improvements |^3| 
and the limit of a large Nc expansion of QCD [Q. The value of the P33 scattering volume was found 
to be smaller than previous determinations. 

6. Off-shell amplitudes for pion-nucleon scattering at low momentum transfers were obtained. The isospin- 
zero combination of the s-wave amplitudes was shown to be relatively large off shell. Rather remarkably, 
the P13, P31 and P33 have the same off-shell dependence for momenta below 300 MeV/c. 

One of us (wrg) gratefully acknowledges discussions with T. E. O. Ericson, B. Loiseau, S. Coon and M. 
Mattis. 
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